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ABSTRACT 

We find that the relative contribution of satellite galaxies accreted at high redshift to the stellar 
population of the Milky Way's smooth halo increases with distance, becoming observable relative to 
the classical smooth halo about 15 kpc from the Galactic center. In particular, we determine line- 
of-sight-averaged [Fe/H] and [a/Fe] in the metal-poor main-sequence turnoff (MPMSTO) population 
along every Sloan Extension for Galactic Understanding and Exploration (SEGUE) spectroscopic line 
of sight. Restricting our sample to those lines of sight along which we do not detect elements of 
cold halo substructure (ECHOS), we compile the largest spectroscopic sample of stars in the smooth 
component of the halo ever observed in situ beyond 10 kpc. We find significant spatial autocorrelation 
in [Fe/H] in the MPMSTO population in the distant half of our sample beyond about 15 kpc from 
the Galactic center. Inside of 15 kpc however, we find no significant spatial autocorrelation in [Fe/H]. 
At the same time, we perform SEGUE-like observations of Y-body simulations of Milky Way analog 
formation. While we find that halos formed entirely by accreted satellite galaxies provide a poor 
match to our observations of the halo within 15 kpc of the Galactic center, we do observe spatial 
autocorrelation in [Fe/H] in the simulations at larger distances. This observation is an example of 
statistical chemical tagging and indicates that spatial autocorrelation in metallicity is a generic feature 
of stellar halos formed from accreted satellite galaxies. 

Keywords: Galaxy: abundances — Galaxy: formation — Galaxy: halo — Galaxy: kinematics and 
dynamics 



1. INTRODUCTION 

The smooth (or classical) halo of the Milky Way 
is the subset of the larger halo population that is 
not in identified substructure that stands out in 
some way from the local background, either photo- 
metrically, kinematically, or chemically. These halo 
substructures have been discovered in star counts, kine- 
matic measure ments, and chemical abun dances (e.g., 
Totten fe Irwin] [l998t iTotten et al.l 12000: llvezic et all 



20001: lYannv et al l 120001: lOdenkirchen et al.l 12001 
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1991 IChiba fc Yoshiil Il99l IHelmi et all Il99l 



Chiba fc Beersl 120001: ikeplev et all 120071: llvezic et all 



2008; IKlement et all 120081: ISeabroke et al.l 1200 
Klement et al.ll2009tlSchlaufman et al.ll2009USmith et al 
Starken burg et al 



20091: llvezic et al 



An et al.1 l2009HSchlaufman et all 120111: iCarollo et all 
20121 ). Despite the large number of substructure de- 



tections, the majority of the stellar content of the halo 
within 20 kpc of the Galactic center is b oth expected 
theor etically (e.g..lMacieiewski et al.ll2011l) an d observed 
(e.g.. lBell et"alll2008l : iSchlaufman et al.ll2009f ) to be in a 
kinematically smooth population. 

While observations of the halo stellar population sug- 
gest that the fraction of the halo in substructure in- 
creases with Galactocentric radius (e.g. , IBell et aTll2008l ; 
Schlaufman et al. 2009; Xuc et al. 2011), the fractions of 
the smooth halo stellar population formed in situ, ac- 
creted in the major mergers that formed the bulk of the 
halo, or accreted in subsequent satellite galaxy minor 
mergers are observationally unconstrained. The ratio of 
in situ star formation to accretion and its dependence 
on distance from the Galactic center may be related to 
the accretion histo r y of Milky Way analog halos (e.g., 
Zolotov et aLll2009l 120101: IHelmi et alll201ll : iFont et all 



20111 ). and the expectation from simulations is that the 



contribution of the accreted com ponen t likel y increases 
with distance (e.g.. IZolotov et al.M2009l 120101: IFont et all 

unni). 

The stellar population in the classical smooth halo 
seems to be kinematically smooth on large scales 
and has three other main properties: (1) it is old, 
(2) its [Fe/H] distribution peaks at [Fe/H] ~ -1.6 
with full width at half maximum > 1 dex, and 
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(3) 
Ca, 



it is 
and 



enhanced in the g-element s O, Mg, Si 



1991 



Ti relative to iron ( e.g.. iRvan fc Norris 



. a Bbl; iMcWilliam et all 119951: lAllende Prieto et al 
20061 ISchlaufman et all 12009ft . iRobertson et al.1 (|2005[ ) 



explained the age, metallicity, a-enhanccment, and kine- 
matic structure of the smooth component of the halo in 
the context of the ACDM model of galaxy formation with 
the accretion of massive M tot ~ 5 x 10 10 Mq halos ~T0 
Gyr in the past. The high mass and short timescale 
for star formation in such massive progenitors of the 
smooth halo are consistent with the observed chemistry. 
In a suite of s i x sim ulations of Milky Way analog halos, 
ICooper et al.l (|2010ft also found that halos acquire the 
bulk of their mass from fewer than fi ve si gnificant con- 
tribut ors. The IRobertson et~aT1 (|2005l ) and lCooper et al.l 
(2010) scenario naturally explains the low [Fe/H] and 
high [a/Fe] in the halo. On the other hand, the halo 
abundance pattern is in contrast to the composition of 
the (surviving) classical dwarf spheroidal (dSph) galax- 
ies, which at the average [ Fe/H] of the smooth halo have 
[a/Fe] closer to solar (e.g.. lMateolll998l;lKirbv et al.ll2010l 
I2011ah . 

In the first paper in this series (jSchlaufman et al.ll200l 
S09 hereafter), we described the results of a system- 
atic, statistical search for elements of kinematically-cold 
halo substructure (ECHOS) in the inner halo metal-poor 
main-sequence turnoff (MPMSTO) population. One by- 
product of the search for ECHOS detailed in IS09I is a 
catalog of MPMSTO stars more than 4 kpc from the 
Galactic plane, more than 10 kpc from the center of the 
Galaxy, within 17.5 kpc of the Sun, and free of both 
surface-brightness and radial-velocity substructure. We 
subsequently refer to this sample as the pure smooth halo 
sample. Though there is likely still very diffuse substruc- 
ture in the catalog, it has been cleaned of substructure to 
the greatest extent possible using existing data. In the 
second paper in this series (|Schlaufman et al.l[20lTL Sll 
hereafter), we analyzed co-added MPMSTO spectra to 
derive the average [Fe/H] and [a/Fe] for ECHOS and for 
the smooth component of the halo along the same line of 
sight as each ECHOS. We found that the MPMSTO stars 
in ECHOS were systematically more metal rich and less 
[a/Fe] enhanced than the MPMSTO stars in the smooth 
component of the halo. We concluded that the chemical 
abundance pattern of ECHOS was best matched by a 
massive dSph galaxy with M to t ^ 10 9 Mq. 

In this third paper of the series, we quantify the de- 
gree of spatial chemical inhomogeneity and spatial varia- 
tion in chemical abundance in the smooth component of 
the halo, using the subs truct ure-cleaned sample of MPM- 
STO stars produced in S09 and the chemical abundance 
technique developed in IS 111 We also investigate the de- 
pendence of these properties on distance from the Galac- 
tic center. 

This concept is analogous to the idea of chemical 
tagging employed in the solar neighborhood. Stars 
formed within the same star forming region have ve ry 
similar abundances (e.g.. iDe Silva et all 12001 l2007al lbl: 
iBubar &: Kingll2"010l : iPo mpcia et alll20ill L and this sim- 
ilarity has led to the use of chemical tagging in con- 
cert with kinematic information to identify groups of 
stars in the solar neighborhood with a common ori- 
gin. The chemical homogeneity of stars formed in 
the same star forming region may be generic (e.g., 



iBland-Hawthorn et al.l 120101 ) . Likewise, many halo sub- 
structures are chemically distinct fro m the kinematically 
smooth halo st ellar population (e.g., Ivezic et all 120081 : 
lAn et al.M2009t iSchlaufman et al.ll2011h . 



To illustrate further, consider the following thought 
experiment. Imagine that a Milky Way analog halo can 
be formed entirely in one of two ways: (1) a scenario in 
which the smooth stellar halo formed in a few dissipative 
major mergers (in which violent relaxation may also have 
been important) and (2) a scenario in which the smooth 
stellar halo formed by the accretion and tidal disruption 
of satellite galaxies (in which diffusion processes are neg- 
ligible). These two limiting cases are vastly simplified 
and cannot accurately reproduce the observed proper- 
ties of the Milky Way. Nevertheless, they are illustrative 
of two opposite, extreme models of halo formation. 

In the stellar halo formed entirely in a few dissipative 
major mergers, the resultant stellar population should be 
very well mixed. Consequently, an observer measuring 
metallicity along many lines of sight closely spaced on 
the sky and in the same distance range throughout the 
halo would find that lines of sight that are closely spaced 
are no more likely to have similar chemical abundances 
than lines of sight that are widely separated. 

In the stellar halo formed entirely by the accretion and 
tidal disruption of satellite galaxies, an observer measur- 
ing metallicity along many lines of sight and in the same 
distance range throughout the halo would find that lines 
of sight that are close together on the sky are also likely 
to be close together in metallicity. This effect arises sim- 
ply because lines of sight that are close together on the 
sky in a small range of distance are more likely to sample 
the stellar debris of only a few recently disrupted satellite 
galaxies. On the other hand, lines of sight that are widely 
separated on the sky would be very unlikely to sample 
the debris of the same accreted galaxy. As a result of this 
effect, the observer would see spatially-correlated metal- 
licities in the halo. 

We measure the degree of spatial chemical inhomo- 
geneity and spatial variation in chemical abundance in 
the smooth component of the halo to differentiate be- 
tween these two possibilities as a function of distance 
from the Galactic center. This paper is organized as fol- 
lows: in Section 2 we define the data we use in this anal- 
ysis. In Section 3, we describe how we determine line- 
of-sight average metallicities in the kinematically smooth 
halo. We also report how we quantify the degree of spa- 
tial autocorrelation in the halo of the Milky Way and 
in halos of simulated Milky Way analogs. In Section 4, 
we discuss the implications of our findings for the forma- 
tion of the Milky Way. We summarize our conclusions in 
Section 5. 

2. DATA 

The Sloan Extension for Galactic Understanding and 
Exploration (SEGUE) survey observed approximately 
240,000 Milky Way stars with apparent magnitudes 
in the range 14 < g < 20.3 with the fiber-fed Sloan 
Digital Sky Survey (SDSS) spectrograph at moderate 
resolution. Spectroscopic targets were selected from 
the combined 11,663 deg 2 ugriz photometric footprint 
of the SDSS and SEGUE. The SDSS telescope and 
spectrograph obtain R 1800 spectra between 3800 
A and 9000 A. The SEGUE data processing pipelines, 
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survey strategy, and radial velocity and atmospheric 
param e ter accur a cies ar e descr i bed in lYannv et aF 



2009), 



2008), 



Smolinski ct al 



paper (jAbazaiian et al.l I200S 
and its instrumentation are 



and the SDSS-II DR7 
The SDSS survey 
described in detail in 



iFukugita et~aH (119961) . IGunn et al.l (11998D. lYork et al.l 



2000 ). jHoggetal 
iPieretal.l (120031). 



2006 |) . iTucker et all 
200. 



(2001) 



ISmithetalJ (120021) . 

IGunn et al 



2004) 



Ivezic et al.l 
2006]), and iPadmanabhan et al 



Way, we perform SEGUE- like observations of both the 
iBullock fc Johnston! (120051) suite of halos and the Via 
Lee et al.l (l2008al(bl. 120111) . lAllende Prieto et al.l Lactea 2 halo (e.g.. [Piemand et al.l 120081: iMadau et al.l 

120081: IZemp et al.l l2009t IRashkov et al J 120111) to deter- 
mine the degree of spatial chemical inhomogeneity and 
spatial variation in chemical abundance present in TV- 
body simulations of Milky Way analog halos. 

3.1. Line of Sight Averaged Chemical Abundances 

For each MPMSTO star in the pure smooth halo sam- 
ple described in Section 2, w e use the M 13 fiducial se- 
quence from lAn et all (|2008f ) to very roughly estimate 
its distance from the center of the Galaxy based on its 
g — r color and apparent r-band magnitude (assuming 
that the Sun is 8 kpc from the center of the Galaxy). 

In order to quantify the spatial variation in chemical 
abundance in the smooth halo (and their dependence on 
distance from the Galactic center) , we subdivide our sam- 
ple into nearby and distant subsamples. The halo stellar 
population may change gradually with distance, and our 
division into nearby and distant regions can only crudely 
capture such a change. However, by splitting the halo 
coarsely into nearby and distant regions, we can compute 
an average value for [Fe/H] and [a/Fe] in two volumes 
along each line of sight. We can then vary the distance 
at which the larger sample is split into the nearby and 
distant subsample in order to explore how the average 
stellar population varies from line of sight to line of sight 
at constant distance from the Galactic center, and how 
the amplitude and spatial scale of that variation change 
with distance. We use three different distance thresh- 
olds to divide our sample into nearby and distant re- 
gions: 12.1, 14.1, and 16.8 kpc. These are the quartilcs 
of the estimated Galactocentric distance distribution of 
our sample. For each line of sight, the result is six over- 
lapping subsamples in distance from the Galactic center. 

We use the spectral co-addition technique summarized 
in Appendix A and described in detail in IS 111 . For each 
subsample and line of sight, we create the "average" 
MPMSTO spectrum for that subsample and line of sight 
by creating a single co-added spectrum from the MPM- 
STO spectra in that subsample and line of sight. Co- 
addition is necessary in this situation because the ap- 
parent faintness of the MPMSTO stars in our distant 
subsamples results in spectral S/N per pixel in individ- 
ual spectra that are not high enough to obtain precise 
and accurate stellar parameter s from the SEGU E Stel- 
lar Parameter Pipeline (SSPP: iLeeet al.N2008allb| , l2fjTTI: 
lAllende Prieto et al.l 120081: iSmolinski et al.l l20lil) . We 
use the [Fe/H] and [a/Fe] reported by the SSPP. In both 
cases, the Fe abundance calculated is an overall metal- 
licity, not just the abundance of iron. 

We verified the accuracy and precision of the co- 
addition procedure using both MPMSTO stars observed 
by SEGUE in the globular clusters M 13 and M 15 and 
SEGUE observations of field MPMSTO stars with stel- 
lar parameters derived from subsequent high-resolution 
spectroscopy. The mean square error (MSE = bias 2 
+ variance) of our co-added abundance analysis ranges 
from 0.1 dex in [Fe/H] at [Fe/H] « -1.0 to 0.20 dex in 
[Fe/H] for [Fe/H] « -1.5. Our MSE estimate includes 
the effect of systematic error, which will not affect the 
spatial autocorrelation analysis we describe in this paper. 
For that reason, these error estimates are conservative 



In IS09L we described the results of a systematic, sta- 
tistical search for ECHOS in the inner halo MPMSTO 
population. We defined the inner halo as the volume 
more than 10 kpc from the Galactic center, within 17.5 
kpc of the Sun, and more than 4 kpc from the Galactic 
plane. While we did find statistically significant evidence 
of substructure along 19 of 137 7 deg 2 lines of sight that 
we searched, on average we found no reason to reject a 
kinematically smooth model for the MPMSTO popula- 
tion in that volume. 

In this paper, we examine the MPMSTO population 
along those 118 7 deg 2 lines of sight where there is no sig- 
nificant surface-brightness substructure and for whi ch we 
found no significant radial- velocity substructure in IS091 
This results in a sample of 9,005 MPMSTO stars in the 
same volume with both photometric and spectroscopic 
data. These MPMSTO stars have both the g - r color 
and significant UV excess expected for the main-sequence 
turnoff of a metal-poor population (for a detaile d de - 
scription of the MPMSTO sample see Section 2 oflS09l). 
Given the magnitude limits of the SEGUE survey, the 
MPMSTO sample was selected because MPMSTO stars 
are the highest-density tracer of the inner halo. At the 
mean g — r color and metallicity of our sample, SEGUE 
radial velocities are precise to 5 km s _1 at g = 18 and 20 
km s _1 at g = 20.3. 

This sample is the largest spectroscopic sample of halo 
stars in the smooth component beyond 10 Galactocentric 
kpc ever assembled. To the extent possible with current 
data, it has been cleaned of both surface-brightness and 
radial-velocity substructure. As such, it allows for an 
unprecedented analysis of the chemical properties of the 
classical smooth halo. We plot in Figure[T]the spatial dis- 
tribution of lines of sight populated by MPMSTO stars 
in the pure smooth halo population. 

3. ANALYSIS 

For each SEGUE line of sight in our pure smooth halo 
sample described in Section 2, we compute the average 
[Fe/H0 and [a/Fe] of the MPMSTO population along 
that line of sight as a function of distance from the 
Galactic center. We subdivide our sample into nearby 
and distant subsamples and compare the average [Fe/H] 
and [a/Fe] in these subsamples to investigate the ef- 
fect of distance from the Galactic center on the chem- 
ical properties of the smooth halo. We quantify the 
inhomogeneity and spatial variation in average chem- 
ical abundances by measuring the spatial autocorrela- 
tion of the line-of-sight-averaged chemical abundances. 
To better understand our observed result for the Milky 

9 [Fe/H] is the context of this paper is the total metallicity as- 
suming all metals are coupled, not truly the iron abundance. 
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upper limits on our errors. 

We give the result of these calculations in Tables [1] 
and [2j We plot the line-of-sight-averaged [Fe/H] dis- 
tribution from Table [1] in Figures [2] and We plot 
the distribution of the data in Tables [T] and [3] in 
the [Fe/H]- [a/Fe] plane in Figure GO Close to the 
Galactic center, the halo is dominated by the classi- 
cal — 1.8 < [ Fe/H] < —1.5 and [a /F el-enhanced pop- 
ulation (e.g.. iRvan &: Norrisl ll991alTbl ; iMcWilliam et al.l 
Il995t lAllende Prieto et al.l l2006i r~The spread in [a/Fe] 
is consistent with a homogeneous [a/Fe] w 0.3 given 
our errors. We find —1.8 < [Fe/H] < —1.5 and 
[a/Fe] « 0.3 in the halo MPMSTO population at 10 
kpc. These values are consistent with high-resolution 
abun dance analyses of loca l halo stars at 100 pc 
'e.g.. lEdvardsson et all 119931: INissen fc Schusterl 119971: 



Hanson et all 119981: iFulbri ght 2000; Prochask a et al 
20001: IStephens fc Boesgaardi 1 2003 IBensbv et al.l 1200 
Reddv et al.ll2003HVenn et alJl2004h . 

Farther from the Galactic center, the halo has a non- 
negligible dSph-li ke component with [Fe/H ] sa —1.0 and 
[a/Fe] « 0.1 fe.g- IKirbv eFalllMol l2011allbh . This com- 
ponent is found preferentially at \b\ < 30°, indicating a 
possible association with the low- latitude substructure 
first identified in Mono c eros (e.g.. iNewberg et alj|2002t 
Belokurov et al.l l2006b| Uuric et al l 120081: llv ezic et al l 
20081: ISchlaufman et al.ll20Tl¥ and perhaps with the d sk 



of the Milky Way. The chemical abundances we find in 
the smooth halo at low Galactic latitude are very similar 
to the chemical abundances we derived in IS11I for this 
low-latitude substructure. This low-latitude substruc- 
ture is not apparent in [a/Fe] alone, as the precision 
of our measurements is much worse in [a/Fe] than for 
[Fe/H]. The observation of this chemical abundance pat- 
tern exclusively at low Galactic latitude indicates that 
it may be related to satellite-induced disk flaring (e.g . , 
IKazantzidis et~alll2008l l200l IPurcell et^l2Wl^01lf) . 

Table[3]provides the mean [Fe/H] of the smooth halo in 
each Galactocentric distance slice, both in the entire sam- 
ple and only for lines of sight with | 6| > 30°. From both 
our current data and our analysis in lSlll . we have reason 
to believe that the known substructure at low Galactic 
latitude is metal rich and quite prominent at the largest 
Galactocentric distances probed by our MPMSTO sam- 
ple. For that reason, we focus on the mean [Fe/H] above 
\b | = 30°. In that case, there are no instances in which 
mean [Fe/H] inferred for the inner subsample differs sig- 
nificantly from the mean [Fe/H] inferred for its comple- 
ment outer subsample. The largest difference we measure 
is for the volume 10 < d < 16.8 kpc and its complement 
volume 16.8 < d < 17.5 kpc: [Fe/H] = -1.77 ± 0.12 
versus [Fe/H] = — 1.53 ± 0.13. That difference is only 
significant at the 1.4er level. As such, it can only be 
regarded as a very tentative hint that there may be a 
positive metallicity gradient in the halo. As a result, we 
find no significant [Fe/H] gradient in the volume probed 
by our MPMSTO sample. 

We have estimated the metallicity bias of our UV ex- 
cess selection for the MPMSTO sample. Relative to 
stars with [Fe/H] = —1.6, stars at [Fe/H] = —3 are 
twice as likely to be selected for spectroscopic follow- 
up, while stars at [Fe/H] = — 1.1 are half as likely to be 
selected. The selection probability follows a smooth, ap- 
proximately parabolic relation over that range in metal- 



licity, and cuts off rapidly below [Fe/H] = — 1. Modeling 
this selection bias and including our measurement errors, 
for a halo population with tr ue mean [Fe/H] = —1 -6 and 
standard deviation 0.4 (e.g.. IRvan fc Norris1ll991a HbT) we 
expect to measure a mean value of [Fc/H] = —1.79, in 
good agreement with the means in Tabic [3] All of the 
metallicity distributions in Figure [3] are consistent with 
a majority population with mean near [Fe/H] = —1.6, 
in agre ement with the measurements of IRvan fc Norria 
(|1991allbT ) plus an additional population of more metal- 
rich stars. The selection bias in our sample is too mild 
for the true population means to be significantly more 
metal rich than about [Fe/H] = —1.5. 

3.2. Simulated Observation of Theoretical Models 
3.2.1. Theoretical Models Used 

We simulate S E GUE- like observations of the 
iBullock fc Johnston] (|2005l ) suite of halos as well as 
the Via Lactea 2 halo. The IBullock fc Johnston! (|2005ft 
simulations used a "hybrid" approach in which they 
modeled accretion events onto a Milky Way-like galaxy 
with A-body particles while modeling the rest of the 
galaxy with a smoothly evolving spherical potential and 
an analytic approximation for dynamical friction. They 
neglected both interactions between particles and the 
gravitational response of the rest of the galaxy to those 
particles. They selected accretion events from merger 
trees for each of their 11 M v i r = 1.4 x 10 12 M galaxies 
in a universe with f2 TO = 0.3, Ha = 0.7, H^h 2 = 0.024, 
h = 0.7, and ag = 0.9. They populated the dark matter 
halos they simulated with star particles using a simple 
star formation law and initialized their orbital planes 
with random orientations. To maximize the computa- 
tional efficiency of their simulations, they only simulated 
accretion events with significant stellar content. This 
approximation resulted in a satellite mass floor of 
Mvh- ~ 10 8 Mp). Given their lack of a s elf-consistent 
potential, the IBullock fc Johnston] (|2005l ) models are 
less able to accurately track the dynamical evolution 
of debris accreted before its host galaxy's last major 
merger. The inner halo is preferentially formed from 
ancient accre tion events, so the p arame tric potentials 
used in the IBullock fc Johnstonl ([2005) models may 
not accurately reproduce the phas e-space properties of 
the in ner halo within 20 kpc. The Bull ock fc Johnston 
(2005) mo dels are further des cribed in IRobertson et al 
( 2005) and lFont et ail (|2006allbh . 

Th e V ia Lactea 2 d ar k mat ter halo (iDiemand et all 
20081 iMadau et all [200l iZemp et all [20091 



Rash kov et al. I 120 111) was simulated using the par- 



allel tree A-body code PKDGRAV2 in a universe with 
Q m = 0.238, Q A = 0.762, h = 0.73, and ct 8 = 0.74. 
The simulation considered the formation of a halo with 
M vir = 1.9 x 10 12 M Q at z = resolved with 1.1 x 10 9 
particles of mass 4, 100 M Q and gravitational softening 
length 40 pc. The 1% most tightly bound dark matter 
particles in each subhalo with M t id > 10 7 M & at infall 
were tagged as star particles in post processing. The 
stellar mass and metallicity of each star particle was 
assigned based on an empirical model that ensures that 
the remnant population of dSph surviving to the present 
day reproduces the obs erved properties of the Milky 
Way's satellite galaxies (Rash kov et all 120111) . The Via 
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Lactea 2 simulation fully accounts for gravitational 
interactions between particles and self-consistently 
determines the growth of the mass of the main halo, 
so it better tracks ancien t accretion events than the 
iBullock fc Johnston! (|2005[ ) models. 

3.2.2. Description of Our Simulated Observations 

To simulate a SE GUE-like observation of one of the 
IBullock fc Johnston! (|2005T ) halos or the Via Lactea 2 
halo, we first place an observer at 8 kpc from the cen- 
ter of the halo along the re-axis of the simulation. This 
choice is arbitrary, as there are no stellar disks in either 
simulation. We consider only those star particles with 
\z\ > 4 kpc that are more than 10 kpc but less than 20 
kpc from the center of the halo. We compute the radial 
velocity and Galactic coordinate of every star particle in 
that volume as it would be viewed by the observer. For 
every SEGUE line of sight listed in Tables Q] and we 
compute the range of I and b observed in the MPMSTO 
star population along that line of sight. We use that 
range of Galactic coordinates to select the star particles 
from the simulations that would have been observed by 
SEGUE along that line of sight. As before, we compute 
the luminosity-weighted mean metallicity of the stars in 
six distances ranges. In that way, we obtain line-of-sight- 
averaged metallicities that match the line of sight aver- 
aged metallicities observed in the Milky Way's smooth 
halo and presented in Tables [T] and [2] One caveat i s that 
the metallicities reported in [Bullock fc Johnston! (j2005l ) 
and iRashkov et al.l ()2011h are calibrated to match the 
true iron abundances of dSph galaxies, and not the total 
metallicities used in this analysis. 

In the right-hand panel of Figure [4] we plot 
the resul t of a SEGUE-lik e obse rvation of halo 
10 from IBullock fc Johnstonl (|2005D . Though the 
IBullock fc Johnston! (|2005l ) mo d els lack the gas-rich ma - 
jor mergers at high redshift that iRobertson et al.l (|2005l ) 
advocated as the origin of the kinematically smooth halo 
stellar population, they do illuminate the properties of a 
stellar halo built entirely from disrupted satellite galax- 
ies. I n an effort to better i s olate the smooth component 
of the IBullock fc Jo hnston (2005) models, we attempted 
to find and remove ECHOS from those models. However, 
there are generally too few star particles to differentiate 
between substructure and smooth component, and our 
search produced no results. For that re a son, w e treat 
an entire halo from IBullock fc Johnstonl ()2005l ) within 
20 kpc of the center of the halo as a smooth compo- 
nent. The star particle tagging prescription applied to 
Via Lactea 2 did not track [a/Fe], so we cannot make 
the equivalent plot for Via L actea 2. 

Figure H] indicates that the IBullock fc Johnstonl (|2005f ) 
models are a poor match to our observation of the smooth 
halo, though there is a hint that they become a better 
ma tch farther from the Galac tic center. In particular, 
the IBullock fc Johnstonl (|2005l ) models have a significant 
[Fe/H] pa —1.0 and [a/Fe] pa 0.1 population in all dis- 
tance bins, while the [Fe/H] pa —1.0 and [a/Fe] pa 0.1 
population we observe in the Milky Way is confined to 
the region far from the Galactic center and at low Galac- 
tic latitude. Indeed, given that our errors in average 
[Fe/H] are at most 0.2 dex, there is no way to recon- 
cile the left hand panels of both plots in Figure [4] with 
measurement error alone. If the low-latitude dSph-like 



composition we observe in the smooth halo at large ra- 
dius and low Galactic latitude is somehow related to the 
stellar disk of the Milky Way, then it is unreasonable to 
expec t it to be reproduced by the IBullock fc Johnstonl 
(2005) models (which have no stellar disk). 

The absence of a significant [Fe/H] pa — 1.0 and [a/Fe] 
pa 0.1 populati on in the smooth hal o of th e Milky Way 
relative to the IBullock fc Johnstonl (|2005h models has 
two main implications. First, that most of the Milky 
Way's halo was accreted at high redshift, and there fore at 
high [a/Fe] as inferred bv IRobertson et all ((2005). Sec- 
ond, that the stellar mass that has phased mixed into 
the smooth halo from massive substructures (with high 
[Fe/H]) is dwarfed by the classical smooth halo compo- 
nent. This observation is in contrast with the chem- 
ical abundances of ECHOS studied in IS 111 where we 
found that more recently accreted stars still kinemati- 
cally distinct as ECHOS have both higher [Fe/H] and 
lower [a/Fe] than the smooth halo. 

Figure [5] shows th e [Fe/H]-[a/Fe] morpholog y of the 
other 10 halos from IBullock fc Johnstonl (120051) . W hile 
each of the 11 halos from Bullock & Johnston (200^) has 
a unique accretion history, the [Fe/H]— [a/Fe] morphol- 
ogy present in Figured] is common to all 11 halos. We 
plot the metallicity distrib utions of both halo 10 from 
IBullock fc Johnstonl ()2005l ) and the Via Lactea 2 halo in 
Figure [6] The strongly-peaked metallicity distribution 
we observe in the inner 20 kpc of the Via Lactea 2 halo 
is the result of two massive satellites accreted at high red- 
shift. These two satellites contributed 80% of the stellar 
mass in that volume and had [Fe/H] = —0.94 and [Fe/H] 
= —1.04. While we showed that we are capa ble o f ac- 
curately measuring [Fe/H] near [Fc/H] = — 1 in lSlll . our 
observational data are inconsistent with a halo having a 
large fraction of its stars at [Fe/H] = — 1 as suggested by 
Figure [6] even after modeling the metallicity bias of the 
UV excess selection. 

3.3. Spatial Autocorrelation 

Tables [T] and [2] give the spatial distribution of [Fe/H] 
and [a/Fe] in the kinematically smooth halo MPMSTO 
population. As we showed in Figure |4j the low-latitude 
substructure at [Fe/H] pa —1.0 and [a/Fe] pa 0.1 is very 
prominent at low Galactic latitude. At the same time, 
Figure 0] does not reveal whether lines of sight that are 
close together on the sky also have similar metallicities. 
As we described in Section 1, a halo that was formed in a 
dissipative process will not have the property that lines of 
sight that are close together on the sky also have similar 
metallicities. On the other hand, a halo that was formed 
from disrupted satellite galaxies will have the property 
that lines of sight that are close together on the sky also 
have similar metallicities. 

To determine whether lines of sight that are close to- 
gether on the sky a lso have simi lar metallicities, we use 
Moran's I statistic (|Moranl ["l950]) to quantify the spatial 
autocorrelation of the chemical abundances: 

n ET E" j fa ~ £) ( X J ~ £) m 

In Equation [T] n is the number of lines of sight, Xi is 
average metallicity inferred for a given line of sight, X 
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is the mean metallicity of the sample, and Wij is the 
spatial weight. In this case, each Wij is the inverse of 
the angular distance between line of sight i and line of 
sight j. The expected value and variance of Moran's / 
under the null hypothesis of no spatial var iation are (e.g., 
ICliff fc Ordlll98lTlWarler fc Gotwavll200l : 



E(7) = 
Vax(J) = 



1 



n-1' 
n 2 S\ — 11S2 



35 2 



(n-l)(n + l)^ 



n- 1 
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(3) 

(4) 
(5) 
(6) 

(7) 

(8) 



Positive values of Moran's / statistic indicate posi- 
tive spatial autocorrelation, while negative values indi- 
cate negative spatial autocorrelation. In other words, if 
lines of sight that are close together on the sky are sig- 
nificantly more likely to have similar metallicities than 
expected based on chance, then Moran's I statistic will 
be significantly positive. 

Though Moran's / statistic gives a good indication of 
the existence of spatial autocorrelation in data, it does 
not indicate the angular scale at which spatial autocor- 
relation is strongest. To constrain the spatial scale of 
autocorrelation, we define the statistic &Ai- 



&AA = 



1 



£ 

jGAj^i 



{Xi Xj ) . 



(9) 



where Xi and the average metallicity inferred for 

the lines of sight indexed by i and j, A is the set of lines 
of sight within a given angular scale of Xi , and m is the 
number of lines of sight in A. 

For each line of sight in Tables [T] and El we compute 
Q A ,i at each of ten angular scales 10°, 20°, ... , 100°. To 
estimate both 6,4 and our uncertainty in &a at each 
angular scale, we use a Monte Carlo simulation. Instead 
of using the set A of all lines of sight within a given 
angular scale of the line of sight denoted by Xi , we create 
a bootstrap set A' by sampling with replacement m lines 
of sight from all m available in A. We compute ®A',i m 
this way 100 times for each line of sight and each angular 
scale. We then compute the mean E(®A',i) and variance 
Va,r(®A',i) of ®A',i for each line of sight and each angular 
scale over the 100 bootstrap realizations of A' . We then 



compute E(@a) by averaging E(®A',i) over all n lines 
of sight and Var(0^) by averaging Va.r(&A',i) over all n 
lines of sight and dividing by y/n. 

To determine the distribution of Qa under the null 
hypothesis that there is no spatial autocorrelation in the 
data, we use a second Monte Carlo simulation. We com- 
pute 6 as before, though instead of using the set A 
of all lines of sight within a given angular scale of the 
line of sight denoted by Xi , we randomly sample m lines 
of sight from all n available regardless of spatial prox- 
imity and call this set B' . We compute ®B',i in this 
way 100 times for each line of sight and each angular 
scale. We then compute the mean E(®B',i) and variance 
Vax(®B',i) of Qb,i for each line of sight and each angular 
scale over the 100 bootstrap realizations of B' . We then 
compute the mean E(0s) and variance Var(Os) under 
the null hypothesis of no spatial autocorrelation by aver- 
aging E(0B',i) over all n lines of sight and by averaging 
Var(6_B',j) over all n lines of sight and dividing by y/n. 

The significance of spatial autocorrelation at a given 
angular scale is then: 



E(e 5 )~E(e j4 ) 

y/Vai(e A ) + Var(9 B ) ' 



(10) 



Like Moran's I statistic, positive values of il indicate spa- 
tial autocorrelation; negative values of £1 indicate spatial 
dispersion. We report Moran's / statistic for the data 
in Table [T] in Tables [4] and [5] we report the importance 
of spatial autocorrelation graphica lly in Figure [71 We 
report Moran's I statistic for the iBullock fc Johnston! 
([20050 halos and Via Lactea 2 in Tables [6] and re- 
spectively; we report the importance of spatial autocor- 
relation graphically in Figure El We find no significant 
spatial autocorrelation in [a /Fe] , likely because our mea- 
surement of that quantity is much less precise than our 
measurement of [Fe/H]. For that reason, we do not con- 
sider [a/Fe] any further. 

We find that there is significant spatial autocorrela- 
tion in [Fe/H] in the Milky Way beyond about 12 kpc 
from the Galactic center, indicated by the very small p- 
values in Table 0] and large positive values of VI in the 
left-hand panel of Figure [7] The spatial autocorrelation 
in [Fe/H] close to the Galactic center is very likely due to 
the prominent substructure at low Galactic latitude with 
[Fe/H] ss —1.0 and [a/Fe] w 0.1. Indeed, if we restrict 
our sample only to lines of sight with Galactic latitude 
I b\ > 30°, the signal of spatial autocorrelation in [Fe/H] 
closer to the Galactic center disappears. This absence is 
indicated by the large p- values in Table [5] and the small 
values of Q, in the right-hand panel of Figure [JJ 

Even at |6| > 30° (i.e., ignoring the low latitude sub- 
structure) , there is still significant spatial autocorrelation 
in [Fe/H] far from the center of the Galaxy in the classical 
smooth halo —1.8 < [Fe/H] < —1.5 and [a/Fe] -enhanced 
population. Given the poor precision of our MPMSTO 
star distance estimator, our constraint on the distance at 
which spatial autocorrelation becomes apparent is weak. 
Nevertheless, the p- values from Moran's / statistic in Ta- 
ble [5] indicate that the effect becomes apparent beyond 
about 14 kpc but before about 16 kpc. The effect is 
very strong at the largest distances probed by or MPM- 
STO star sample. The smallest angular scale that we 
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are sensitive to spatial autocorrelation in [Fe/H] is 10°. 
Interestingly, in Figure [7] the statistic il peaks between 
20° and 30°. This is the characteristic angular scale of 
spatial autocorrelation. At 14 kpc - the median Galac- 
tocentric distance of the MPMSTO stars in our sample 
- that corresponds to a physical scale of approximately 5 
kpc. This is much larger than the scale of the disrupted 
satellite galaxies that contributed the stellar populations 
that are likely the source of this effect. 

After we identified the characteristic angular scale of 
spatial autocorrelation in metallicity, we went back to the 
data in Table [T] and computed Moran's / statistic in 30° 
windows around each line of sight. While we found that 
a handful of lines of sight displayed significant spatial au- 
tocorrelation in their own local 30° window, the p- values 
that resulted from those observations were orders of mag- 
nitudes larger than the p- values we computed for the full 
sample of 118 lines of sight. Many lines of sight display 
no significant spatial autocorrelation in their own 30° 
window. Consequently, the signal of spatial autocorrela- 
tion that we observe is not the result of a small number 
of discrete locations in the sky. Instead, the significant 
spatial autocorrelation in [Fe/H] we observe is coming 
from everywhere in the Galactic halo. 

We also recomputed Moran's / statistic using the sig- 
nal to noise-weighted line-of-sight average metallicities 
created by averaging DR8 SSPP metallicities for all in- 
dividual stars in the same lines of sight and volumes 
described in Table [1] That analysis produced a quali- 
tatively similar result to that described above, albeit at 
much lower statistical significance. 

We do not identify any significant spatial autocorre- 
lation in [Fe/H] in the M ilky Way analog halos from 
iBullock fc Johnston! ([20051 ) in the volume between 10 and 
20 kpc from the halo center. This lack of spatial auto- 
correlation is indicated by the large p- values in Table [6] 
and small values of ft in the left-hand panel of Figure [SJ 
The large p- values and small f2 values we present for 
halo 10 from IBullock fc Johnston! (|2005[ ) are representa- 
tive of the values we obtain from eq uivalent observations 
of the other ten halos presented in IBullock fc Johnston! 
(2005). The smoothness of the star-particle distribution 
near the center of the halo could be physical, or it could 
be a byproduct of the limitations of the simulations. 
We do observe significant spatial autocorrelation in the 
IBullock fc Johnston! (|2005l ) halos beyond 20 Gala c tocen - 
tric kpc though, where the IBullock fc Johnston! (2005) 
simulations likely model halo formation more accurately. 

On the other hand, we do observe significant spatial 
autocorrelation in [Fe/H] in the Via Lactea 2 halo in the 
volume between 10 and 20 kpc from the halo center. The 
significant spatial autocorrelation in [Fe/H] is indicated 
by the small p- values in Table [7] and large values of ft 
in the right-hand panel of Figure [8j Since we observe 
significant [Fe/H] spatial autocorrelation in the inner re- 
gions of the Via Lactea 2 halo, but little [Fe/H] spatial 
auto correlation in the inn e r regi ons of any of the 1 1 halos 
from[Bullock fc Johnston! ([20051 ). the Via Lactea 2 obser- 
vation is not likely due to it s uniqu e accretion history. 

The IBullock fc Johnston! (|2005l ) halos and the Via 
Lactea 2 halo have similar mass and were simulated 
using similar cosmological parameters. One small dif- 
fer ence between the two simu lations is the value of as, 
as IBullock fc Johnston! (|2005l ) used cr 8 = 0.9 while Via 



Lactea 2 used eg = 0.74. While this affects the abun- 
dance of Milky Way analog halos, it should not much af- 
fect t he properties of i ndividual Milky Way analog halos 
(e.g.. lFont et al.ll2011[ ). As a result, the two simulations 
should not have statistically disparate accretion histories. 
For that reason, the fact that we do not observe signifi- 
cant spatial autocorrelation in [Fe/ HJ in the inner regions 
of any of the 11 halos presented in IBullock fc Johnston! 
(2005) , but do observe the effect in the inner region of the 
Via Lactea 2 halo, indicates that the discrepancy is not 
likely to result from differences in halo merger history. 

The observed differences between the two calculations 
may be a result of the higher resolution and cosmolog- 
ically self-consistent potential of the more modern Via 
Lactea 2 simulation. Another hypothesis is that ac- 
cretion along a filament in the Via Lactea 2 halo pro- 
duced the observed spatial autocorr elation, as the orbital 
planes of infalling satellites in the Bul lock fc Johnston! 
(2005) models were chosen randomly and independently. 
The simulations necessary to verify the latter speculative 
hypothesis would be very computational intensive, and 
therefore are beyond the scope of this analysis. 

4. DISCUSSION 

We have three observations that a complete theoreti- 
cal model for the formation of the inner halo must re- 
produce: (1) the halo metallicity distribution, (2) the 
joint halo [Fe/H]-[a/Fe] distribution, and (3) the signif- 
ican ce of spatial autocorrelatio n in [Fe/H] in the halo. 
The IBullock fc Johnston! (|2005[ ) models to not reproduce 
any of the three measurements, as expected by those au- 
thors because of their simplifying assumptions. The Via 
Lactea 2 model does match (3) even as it fails to repro- 
duce (1) and (2). While self-consistent, high-resolution 
pure-accretion models may be able reproduce (3), they 
have not yet been shown capable of reproducing (1) and 

As described in Section 3.3, we do not observe signif- 
icant spatial coherence in [Fe/H] in the halo of Milky 
Way inside of 15 kpc. This observation is consistent 
with both a smooth halo formed through a combination 
of in situ star formation and dissipative major mergers 
at high redshift and a smooth halo formed through the 
phase mixing of the debris of disrupted s atellite galaxies. 
However, the IBullock fc Johnston! ([200 51) halos also sim- 
ulate the chemical abundance structure of a halo built 
entirely from accreted galaxies. As we observed in Fig- 
ures mm the overall metallicity of such halos are incom- 
patible with the overall metallicity of the halo of the 
Milky Way. As a result, our observation of the chemi- 
cal structure and spatial autocorrelation in [Fe/H] in the 
smooth halo supports the idea of in situ star formation 
and dissipative major mergers at high redshift as the 
origin of the smooth ha l o (e.g.. iRobertson et all 120051; 
Zolotov et al.|[200l [20101: iFont et al.l l2011l:lTissera et al.l 
2012tlMcCarthv et al.ll201ll) ~ 

We observe significant spatial autocorrelation in [Fe/H] 
in the halo of Milky Way exterior to about 15 Galacto- 
centric kpc. The spatial autocorrelation is a global prop- 
erty of the entire halo beyond 15 kpc, and it cannot be 
simply associated with a few locations on the sky. These 
observations favor a halo model in which the relative con- 
tribution of accreted galaxies to the stellar population of 
the halo increases with radius, becoming observable in 
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our data relative to the smooth halo population at about 
15 kpc. 

The dominance of the classical smooth halo close to 
the center of the Galaxy does not mean that disrupted 
satellite galaxies have not contributed to the halo stellar 
population at small distance. The most massive galaxies 
accreted by the Milky Way will quickly spiral into the 
center of the halo due to dynamical friction and deposit 
most of their stars in the inner halo. The chemical abun- 
dances of the recently accreted stars in the inner halo in- 
dicate that they were likely formed in the most massive of 
the s atellite galaxies recently accreted by the M ilky Way 
(e.g-. lFont etaLll2006bl : ISchlaufman et al.ll201lD . Conse- 
quently, there is likely a substantial number of accreted 
stars close to the Galactic center. However, the classical 
smooth component of the halo peaks close to the Galactic 
center too, dominating the component of the halo from 
tidally-disrupted satellite galaxies. 

Though the densities of both the classical smooth halo 
and its tidally-disrupted satellite galaxy counterpart de- 
crease with distance, our observation that there is signifi- 
cant spatial autocorrelation in [Fe/H] far from the center 
of the Galaxy indicates that the classical smooth com- 
ponent decreases more quickly with distance than the 
tidally-disrupted satellite galaxy component. Our con- 
clusion is that 15 kpc from the center of the Galaxy is 
the approximate radius at which the accreted halo be- 
comes non-negligible in comparison to the halo compo- 
nent formed in a combination of in situ star formation 
and dissipative major mergers at high redshift. While 
our observation of spatial autocorrelation in [Fe/H] in- 
dicates that the contribution of accretion to the smooth 
component becomes important at about 15 kpc from the 
Galactic center, it is difficult to quantify. In IS09I . we 
found that about 30% of the inner halo volume has 10% 
of its MPMSTO population in diffuse ECHOS. We spec- 
ulate that this diffuse ECHOS population may be the 
origin of the signal we observe. 

Our observation of spatial autocorrelation in [Fe/H] is 
a first example of in situ statistical chemical tagging in 
the halo. Our result suggests that future surveys of the 
halo stellar population using high-resolution multi-object 
spectroscopy will identify dynamically ancient groups of 
stars with similar abundance patterns indicative of for- 
mation in a single star forming region. In particular, 
future surveys with the ability to chemically tag stars in 
situ more than 15 kpc from the Galactic center should 
find many prominent chemical substructures with angu- 
lar scales comparable to 30°. 

The observation that the classical smooth halo compo- 
nent ceases to be the dominant component of the stellar 
population of the halo at about 15 kpc from the Galac- 
tic center is consistent w i th the dua l-halo hypothesis 
advanced in iCarollo et"al~l (|2007l l2010f > and iBeers et all 
(12011) . Though we do not observe a significant metallic- 
ity gradient in the smooth halo, our observations do not 
preclude the existence of a metall icity gradient. Like- 
wise, this result is analogous to the lde Jong et al.l (|2010f ) 
observation that the halo density profile falls off rapidly 
with distance from the Galactic center to a distance of 
about 15 to 20 kpc. Beyond that distance, they found a 
substa ntially lower densit y, slowly varying halo density 
profile. ISesar et al.l (|2011f) observed something similar at 
even larger radius. 



5. CONCLUSION 

We find significant spatial coherence in [Fe/H] in the 
MPMSTO population in the stellar halo of the Milky 
Way beyond about 15 kpc from the center of the Galaxy. 
That spatial coherence suggests that the relative contri- 
bution of disrupted satellite galaxies to the stellar popu- 
lation of the smooth halo increases with radius, becoming 
observable in our data relative to the classical kinemati- 
cally smooth halo at about 15 kpc from the Galactic cen- 
ter. SEGUE-like observations of the Via Lactea 2 halo 
indicate that spatial autocorrelation in [Fc/H] is a generic 
feature of stellar halos formed entirely by accretion. We 
find that this spatial autocorrelation is strongest on an- 
gular scales between 20° and 30°, corresponding to a 
physical scale of about 5 kpc at 15 kpc. Though we find 
that a significant fraction of the stellar halo of the Milky 
Way beyond about 15 Galactocentric kpc is likely com- 
prised of the phase-mixed debris of satellite galaxies, the 
morphology of the halo in the [Fe/H]-[a/Fe] plane inside 
of 15 kpc is not well matched by phase-mixed tidal de- 
bris. Rather, the smooth halo inside of 15 kpc is likely 
formed through a combination of in situ star formation 
and dissipative major mergers at high redshift. 
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APPENDIX 

A. DESCRIPTION OF CO-ADDITION ALGORITHM 



Using the co-addition technique described in detail in IS 111 for each subsample we create the "average" MPMSTO 
spectrum for that subsample by creating a single co-added spectrum from all of the MPMSTO spectra in that sub- 
sample. We include in the co-added spectra only those spectra that correspond to MPMSTO stars within the 500 K 
range in effective temperature T e g that maximizes the S/N noise in the resultant co-added spectrum. We shift each 
MPMSTO spectrum eligible for inclusion in a co-add to a heliocentric radial velocity v r — km s _1 . We then use 
natural cubic spline interpolation to interpolate both the spectrum and its inverse variance on to a common grid in 
wavelength. Next, we numerically integrate the area under the curve defined by the spectrum and normalize both 
the spectrum (by dividing by the normalization factor) and the inverse variance (by multiplying by the normalization 
factor squared) to ensure that each spectrum that is to be included in the co-add has the same scale. For each popu- 
lation of interest, we then create an ensemble of realizations of the co-added spectrum by bootstrap resampling from 
the set of radial- velocity zeroed, interpolated, and normalized spectra that belong to that population. Each spectrum 
contributes to each wavelength bin in proportion to its in verse variance in that bin relative to the other spectra se- 
lected for co-addition. We tested the co-addition process in IS 111 and obtained good agreement between globular cluster 
metallicities produced by co-addition and their known metallicities from high-resolution spectroscopy. 

We also determine an equivalent ugriz photometric measurement for our co-add spectra by computing a weighted 
average of the ugriz photometric measurement of the individual stars in each bootstrap co-add, using the mean S/N 
between 3950 A and 6000 A as the weight. 

We use the SEGU E Stellar Parameter Pipeline (SSPP; ILee et all l2008allbl . I20TTI: lAllende Prieto et al.1 120081: 
iSmolinski et al.|[2011l ) to determine [Fe/H] and [a/Fe] from the co-added "average" spectrum. The SSPP uses Sloan 
spectroscopy and ugriz photometry to infer the stellar atmosphere parameters (T g, logg, [Fe/H], and [a/Fe]) of stars 
observed in the course of SDSS and SEGUE. The SSPP implements a multimethod algorithm in which many different 
techniques are used to compute the stellar parameters. The SSPP then averages the result of each method known 
to be valid in a given color and S/N range to determine the final T c g, \ogg, [Fe/H], and [a/Fe] reported for all stars 
observed in the SDSS and SEGUE surveys. 
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Figure 1. ECHO-free pure smooth halo sample line-of-sight spatial distribution in Galactic coordinates. 
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Figure 2. Density of line-of-sight-avcraged metallicity distribution between 10 and 17.5 kpc from the Galactic center in the smooth 
component of the halo kernel smoothed to 0.1 dex. The orange curve is the distribution for all lines of sight, while the blue curve is the 
distribution for lines of sight with \b\ > 30°. The mean [Fe/H] of each subsample is indicated by the vertical dashed line. Our estimates 
of the mean metallicity include both random (~0.01 dex) and systematic error (~0.1 dex). The distribution for the full sample has a 
heavy tail to high metallicity that is not present at high Galactic latitude. As a result, the mean of the high Galactic latitude sample is 
systematically m ore metal poor by about 0.1 dex. This offset is likely related to the metal-rich low-latitude substructure first identified in 
Monoceros (e.g.. INewbere et ah 2002: Bclokurov ct al. 2006b; Juric et alj|200a Hvezic et alj |2008; Schlaufm an et al.ll20Tll 1. 
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Figure 3. Density of line-of-sight-averaged metallicity distribution as a function of distance from the Galactic center in the smooth 
component of the halo kernel smoothed to 0.1 dex. The y-axis scale is the same for each panel, and its numerical range is arbitrary such 
that the area under each curve integrates to 1. We split our tracer sample in two very roughly in distance at three different split points: 
12.1, 14.1 and 16.8 kpc. Each row corresponds to the result for a single split-point. In all cases, the nearby part of the sample is in the left 
column and the distant part of the sample is in the right column. The mean [Fe/H] of each subsamplc is indicated by the vertical dashed 
line. Our mean metallicities are in the range —1.8 < [Fe/H] < —1.5. We give the average [Fe/H] and uncertainties in Table[3] We observe 
no significant metallicity gradient in our sample, though there may be a hint of a positive metallicity gradient at large radius. 
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Figure 4. Line of sight averaged abundance distribution as a function of distance from the Galactic center in the smooth component 
of the halo in the [Fc/H]— [a/Fe] plane. We split our tracer sample in half very roughly in distance at three different split points: 12.1, 
14.1 and 16.8 kpc. Each row corresponds to the result for a single split-point. In all cases, the nearby half of the sample is in the left 
column and the distant half of the sample is in the right column. Left: from Tables [T1 and l2l the result for our pure smooth halo MPMSTO 
population. The typical precision of our measurement is given by the gray error bars. Our precision is a function of metallicity, so if there 
are two sets of error bars present, the left bars apply for [Fe/H] < —1.0 while the right bars apply for [Fe/H] > —1.0. We plot in blue 
lines of sight with |b| > 30° and in orange lines of sight with |b| < 30°. The smooth halo is dominated by the familiar —1.8 < [Fe/H] 
< —1.5 and [a/Fe]-enhanced population, though at larger distances and low Galactic latitude the pure smooth halo also has a dSph-like 
component with [Fe/H] ~ —1. and [a/Fe] m 0.1. The feature at low Galact i c latitude is likely related to the low- l atitude substructure firs t 
identified in Monoceros (e.g.. INewberg et al.ir2003 : IBelokurov et al.M2006bl : lJuric et al.ll200Sl; llvezic et aljIlOOSl; IBchlaufman et al.|[2(7Tll) . 
Right: the equivalen t lum inosity- weighted average abundance distribution that results from a SEGUE-like observation of halo 10 from 
Bullock & Johnston (2005). The result is similar for all 11 halo s presented in that paper. We discuss the significance of the observed 
differences between our observations of the Milky Way and of the Bullock & Johnston (2005) models in Section 3.2. 
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Figure 5. Line of sight averaged abundance distribution in the [Fe/H]-[a/Fe] plane for the other 10 halos from Bullock & Johnston 
(2005). We include only those star particles with \z\ > 4 kpc that are more than 10 kpc but less than 20 kpc from the center of the halo. 



16 



SCHLAUFMAN ET AL. 



[Fe/H] 
-2 -1 



[Fe/H] 
-2 -1 



[Fe/H] 
-2 -1 




-BJ05: 10 kpc < R < 14.1 kpc 



BJ05: 10 kpc < R < 12.1 kpc 



- -BJ05: 14.1 kpc < R < 20 kpc 



LBJ05: 10 kpc < R < 16.8 kpc J LBJ05: 16.8 kpc < R < 20 kpc 




BJ05: 12.1 kpc < R < 20 kpc 










[Fe/H] 



[Fe/H] 1 




[Fe/H] 



[Fe/H] 



Figure 6. Density of linc-of-sight-avcrag ed metalli c ity dis tribution as a function of distance fro m the center of the halo kernel smoothed 
to 0.1 dex. Left: halo 10 from lBullock fc Johnston] H 20051 ). Right: the Via Lactea 2 halo (e.g., I >:< inai.d al. 2008; Mada u~et""aTl 120081 : 
Zcmp et al. 2009). Dark matter particles in the Via Lactea 2 halo were tagged as star p articles with metallicities determined by a simple 
satellite instantaneous enrichment at accretion model described in lRashkov et al. (2011). The strongly peaked metallicity distribution of 
the Via Lactea 2 halo is the result of the accretion of two massive satellites with [Fe/H] m —1.0 at high redshift. Neither dark matter only 
JV-body simulation is a good match for our observations of the metallicity of the smooth halo, indicating the importance of gas and in situ 
star formation in the origin of the smooth halo within 17.5 kpc of the Galactic center. 
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Figure 7. Quantification of spatial autocorrelation in [Fe/H] in the Galaxy's pure smooth halo MPMSTO population as a function of 
Galactocentric distance and angular scale. Positive values of the statistic f! indicate spatial autocorrelation. Left: the result for all lines of 
sight. Right: the result for lines of sight with |b| > 30°. Ignoring the low-latitude substructure and looking only at high Galactic latitude, 
there is no spatial autocorrelation in [Fe/H] interior to 14.1 kpc. Beyond 14.1 kpc however, there is significant spatial autocorrelation in 
[Fe/H] with characteristic angular scale between 20° and 30°. At 14 kpc - the median Galactocentric distance of the MPMSTO stars in our 
sample - that corresponds to a physical scale of approximately 5 kpc. This observation suggests that the relative contribution of disrupted 
satellite galaxies to the stellar population of the smooth halo increases with distance, becoming comparable to the classical smooth halo 
component at about 15 kpc. 
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Figur e 8. Quantification of spatial autocorrelation in [Fe/H] in simulations of Milky Way analog formation. Left: halo 10 from 
IBullock & Johnston (2005). Right: the Via Lactea 2 halo. There is no significant spatial autocorrelation in [Fe/H] within 20 kpc of 
the center of halo 10 from Bullock & Johnston (2005). There does appear to be significant spatial autocorrelation in [Fe/H] further from 
the center of the halo, where the Bullock & Johnston (2005) simulations likely model halo formation more accurately. On the other hand, 
there is significant spatial autocorrelation in [Fe/H] within 20 kpc of the center of the Via Lactea 2 halo. The appearance of significant 
spatial autocorrelation in [Fe/H] with 20 kpc of the halo center and the angular scale of that spatial autocorrelation in the Via Lactea 
2 halo are very similar to the same quantities measured in the halo of the Milky Way. The rapid drop-off in the significance of spatial 
autocorrelation in [Fe/H] at large angular scale in the Via Lactea 2 halo may indicate that there is less large-scale diffuse substructure in 
the Via Lactea 2 halo than in the Milky Way. The extra large-scale substructure in the Milky Way may be related to the ongoing merger 
with the Sagittarius dwarf galaxy. 
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SCHLAUFMAN ET AL. 

Table 1 — Continued 
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a In/out split at d ~ 12.1 kpc 
b In/out split at d ~ 14.1 kpc 
c In/out split at d ~ 16.8 kpc 



THE FORMATION OF THE MILKY WAY THROUGH COLD HALO SUBSTRUCTURE. III. 21 



Table 2 

Spatial and Radial Distribution of [o/Fe] in the Smooth Component of the 

Halo 
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O OI 
U.Zl 


O 1 ft 
U.lo 


o oo 
u.zu 


O OQ 
U.ZO 


OIK 
U.IO 


O 1 T 
U.l ( 


O 1 1 
U.ll 


o 07 

U.Z i 


OIO 
U. 1U 


O 1A 

U.04 


O Ofl 

u.uo 


O OQ 

u.uy 


OIO 
U. 1Z 


O 40 


O Ofl 

u.uo 


o oq 

U.ZO 


O OO 

u.zz 


O 11 
U.OO 


O 00 

u.zz 


oo 
u.zu 


O 00 

u.zz 


O OPx 
U.ZO 


O 00 

u.zz 


o 02 
u.zo 


n iq 
U.lo 


O QO 
U.oZ 


o r\T 

U.U / 


O 07 
U.U / 


OIK 
U.IO 


O QO 
U.oZ 


O 07 
U.U / 


n iq 

U. lo 


OIO 
U. 1Z 


O 11 
U.OO 


O OQ 

u.uy 


oo 
u.uz 


U.IO 


O Q7 
U.O ( 


O Oft 
U.Uo 


n 17 
U. 1 1 


O np 

u.uo 


O Qfl 
U.OO 


O 07 
U.U / 


O Ofl 

u.uo 


OIO 
U. 1U 


O OQ 
U.ZO 


O 07 
U.U / 


o oi 

U.Zl 


OIO 

u.iy 


O Ofl 
U.ZO 


O (17 
U.U / 


O OO 

u.uy 


o oo 
u.zz 


O OO 
U.ZU 


O Ofl 

u.uo 


o qo 
U.oU 


O nP 

u.uo 


O AA 
U.44 


O Ofl 

u.uo 


O 1 A 
U. 14 


OIO 
U.1U 


O QO 

u.oy 


O Ofl 

u.uo 


o i x. 

U. 10 


o oo 
u.uy 


O 1 O 
U. 1Z 


O 1 1 
U.ll 


O Ofl 
U.ZO 


O 1 1 
U.ll 


O O/l 
U.U4 


O 1 o 
U. 1U 


o o^ 
u.uo 


OIO 
U. 1U 






O 0^ 

u.uo 


OIO 
U. 1U 


17 
U.l ( 


17 
U.l* 


U. lo 


o oo 
u.zz 






O 1 Q 

U.lo 


o oo 
u.zz 


o np 
u.uo 


o oo 
u.zz 


U. lo 


o oo 
u.zz 






O 1 ft 
U.lo 


o oo 
u.zz 


O OP 
U.ZO 


O 1 Q 

u. iy 


n it 
u.o / 


O Oft 
U.Uo 


O 1 1 
U.lo 


U. 10 


O QQ 

u.oy 


o oo 
u.zu 


O QS 
U.oo 


OIO 
U. 1Z 


o 07 

U.Z ( 


OIO 
U.1U 


O 1*^ 
U.oO 


O Ofl 

u.uo 


O 1A 
U.o4 


O 1 Q 

U.lo 


O QO 

u.oy 


O Ofl 

u.uo 


o ofi 
U.ZO 


O 07 
U.U i 


O 1A 
U.o4 


o r\T 

U.U / 


O QO 
U.oZ 


OIO 
U.1U 


O QO 

U.oU 


o nT 

U.U / 


o o^ 
u.zo 


OIO 
U. 1U 


O Qfl 
U.OO 


O Ofl 

u.uo 


O OO 

u.uu 


oo 
u.zu 


O IP 

u.oo 


O Ofl 

u.uo 


o in 
U.oU 


o of; 
u.uo 


O 1 7 
U.l* 


O Oft 
U.Uo 


O 1 1 
U.ll 


o oo 
u.uy 


O 1A 
U.o4 


O Ofl 

u.uo 


o 1 a 

U. 14 


O 1 1 
U.ll 


O AA 
U.44 


o oo 
u.uy 


O 1 Q 
U.lo 


O 1 A 
U.14 


O IT 
U.O ( 


O Oft 
U.Uo 


O OS 

u.zo 


O 07 
U.U i 


O AA 


O 07 
U.U / 


O Q4 


O Oft 
U.Uo 


O 41 


O Ofl 

u.uo 


o qo 
u.oz 


O 07 
U.U i 


O 41 
U.^l 


O 07 
U.U / 


O QQ 
U.OO 


O 1 Q 
U.IO 


O Qft 
U.OO 


O 07 
U.U / 


o qi 
U.ol 


o oo 
u.uy 


O QO 
U.oZ 


O Oft 
U.Uo 


O Oft 
U.Zo 


OIO 
U. 1Z 


O AO 
\).£±Z 


O Oft 
U.Uo 


O 07 
U.Z ( 


OIO 
U. 1U 


O 01 
U.Zl 


OIO 
U. 1Z 


O 1 Q 
U. IO 


14 
U. 14± 


O 04 


11 
U.ll 


n 19 


O 00 

u.zz 






OIO 
U. 1Z 


O 00 

u.zz 


14 
U. 1^ 


O 00 

u.zz 


n 1 c; 
U. 10 


O 1 ft 
U.lo 


O /IO 
U.4U 


O 1 o 
U. 1U 


n it; 
U. 10 


O 1 ft 

U.lo 


O 1 ft 
U.lo 


O 1 Q 

u. iy 


o a i 

U.41 


o np 

u.uo 


O QO 

u.oy 


O Ofl 

u.uo 


o An 

U.4U 


O Oft 
U.Uo 


O IT 
U.O ( 


o nT 

U.U / 


o Qft 

U.OO 


14 
u. ±*± 


O Q7 
U.O / 


17 
U.ll 


O 4Q 


14 
U. 1^ 


O QO 
U.OZ 


11 
U.ll 


o ao 

U.'-iZ 


O 07 
U.U i 


O 07 
U.Z/ 


O 07 
U.U / 


O OQ 

u.zy 


11 
U.ll 


O Q^ 
U.OO 


O Ofl 

u.uo 


0.29 


0.07 


0.33 


0.07 


0.21 


0.08 


0.26 


0.07 


0.16 


0.08 


0.35 


0.06 


0.13 


0.18 


0.43 


0.06 


0.35 


0.07 


0.43 


0.06 


0.36 


0.09 


0.32 


0.06 


0.28 


0.06 


0.25 


0.06 


0.27 


0.08 


0.33 


0.06 


0.27 


0.19 


0.19 


0.22 


0.41 


0.12 


0.16 


0.21 


0.36 


0.08 


0.31 


0.06 


0.26 


0.12 


0.44 


0.06 


0.02 


0.21 






0.02 


0.21 






0.24 


0.15 






0.24 


0.15 


0.18 


0.16 


0.43 


0.06 


0.30 


0.06 


0.36 


0.18 


0.35 


0.06 


0.30 


0.06 


0.32 


0.07 


0.38 


0.07 


0.31 


0.07 


0.18 


0.07 


0.30 


0.06 


0.12 


0.11 


0.28 


0.06 


0.33 


0.08 


0.25 


0.07 


0.27 


0.16 


0.23 


0.07 


0.32 


0.15 


0.28 


0.07 


0.17 


0.20 


0.29 


0.07 


0.20 


0.15 


0.38 


0.07 


0.21 


0.17 


0.35 


0.09 


0.29 


0.06 


0.22 


0.06 


0.34 


0.07 


0.19 


0.06 
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Table 2 — Continued 



RA 


Dec 


] 


b 


[JjJldLrL- 


IpleitrG 


frv/Fpl- a 
[uy rej m 


err 




("Hptr"! 
V aL &! 


(deg) 


CHptr~l 

V aL &! 










45.24 


5.74 


171.39 


-44.61 


2307 


2322 








158.57 


44.34 


171.74 


57.63 


2557 


2567 


0.27 





07 


48.25 


5.48 


174.65 


-42.75 


2335 


2340 








51.24 


5.20 


177.71 


-40.80 


2334 


2339 








57.20 


10.31 


178.00 


-33.00 


2679 


2697 








113.49 


40.89 


178.00 


25.00 


2683 


2701 








166.97 


38.59 


178.45 


65.54 


2856 


2871 


0.43 





07 


37.40 


-8.47 


178.72 


-60.22 


2047 


2067 


0.36 





07 


47.00 


0.00 


179.01 


-47.44 


2048 


2068 








111.29 


37.62 


180.89 


22.44 


2941 


2946 








112.51 


35.99 


182.90 


22.87 


2053 


2073 








53.00 


0.00 


184.53 


-42.87 


2049 


2069 








134.44 


37.13 


185.88 


40.31 


2380 


2400 








110.74 


31.44 


187.00 


20.00 


2677 


2695 








152.52 


35.29 


189.36 


54.80 


2387 


2407 


0.26 





Id 


55.43 


-6.41 


193.71 


-44.60 


2050 


2070 








59.42 


-5.86 


195.91 


-40.94 


2051 


2071 








144.00 


30.05 


197.01 


47.32 


2889 


2914 








118.00 


23.20 


197.73 


23.16 


2891 


2916 








116.00 


18.18 


201.97 


19.54 


2054 


2074 








71.40 


-5.68 


203.00 


-30.48 




2942 








152.38 


25.93 


205.39 


53.92 


2386 


2406 


0.22 





15 


156.53 


17.74 


220.87 


55.27 


2853 


2868 


0.40 





06 


141.56 


7.30 


225.30 


37.58 


2382 


2402 








169.09 


19.29 


227.63 


66.83 


2857 


2872 


0.32 





07 


127.96 


-4.33 


229.00 


20.00 


2807 


2828 








156.63 


8.82 


234.18 


51.20 


2854 


2869 


0.25 





07 


150.00 


0.00 


239.10 


40.72 


2852 


2867 


0.13 





21 


181.81 


19.97 


245.85 


77.61 


2893 


2918 


0.36 





06 


168.77 


9.61 


245.98 


61.30 


2393 


2413 


0.13 





07 


162.00 


0.00 


250.28 


49.82 


2389 


2409 


0.40 





07 


174.00 


0.00 


266.09 


57.37 


2862 


2877 


0.20 





07 


180.94 


9.98 


267.43 


69.50 


2963 


2965 


0.20 





06 


167.16 


-16.21 


270.00 


40.00 


2690 


2708 


0.37 





12 


169.73 


-11.87 


270.00 


45.00 


2859 


2874 


0.28 





09 


172.22 


-7.52 


270.00 


50.00 


2861 


2876 


0.37 





07 


181.00 


0.00 


278.20 


60.57 


2892 


2917 


0.35 





07 


186.00 


0.00 


288.15 


62.08 


2558 


2568 


0.16 





11 


189.00 


0.00 


294.52 


62.62 


2895 


2920 


0.26 





08 


191.00 


-2.50 


299.18 


60.32 


2897 


2922 


0.13 





10 


191.16 


-7.83 


300.00 


55.00 


2689 


2707 


0.18 





09 


198.00 


0.00 


314.09 


62.43 


2901 


2926 


0.22 





14 


194.57 


19.74 


315.26 


82.45 


2899 


2924 


0.25 





06 


205.28 


9.39 


338.75 


68.73 


2903 


2928 


0.15 





11 


217.40 


8.47 


358.72 


60.22 


2908 


2933 


0.30 
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/FpI * a 


err 


[Uij rcj ln 


err 


\n /FpI t b 
[<-<•/ L "J out 


err 


In /FpI • c 


err 


0.17 


0.08 


0.30 


0.11 


0.28 


0.20 


0.27 


0.10 


0.28 


0.06 


0.22 


0.07 


0.24 


0.06 


0.28 


0.06 


0.29 


0.16 


0.18 


0.21 


0.38 


0.21 


0.33 


0.15 


0.15 


0.19 


0.52 


0.16 


0.11 


0.18 


0.13 


0.20 


0.33 


0.22 






0.33 


0.22 


0.37 


0.22 


0.16 


0.22 






0.16 


0.22 






0.39 


0.06 


0.34 


0.07 


0.45 


0.08 


0.42 


0.06 


0.21 


0.08 


0.27 


0.07 


0.16 


0.11 


0.28 


0.07 


0.35 


0.08 


0.35 


0.07 


0.29 


0.15 


0.34 


0.07 


0.21 


0.20 






0.21 


0.20 






0.15 


0.14 






0.15 


0.14 






0.38 


0.07 


0.26 


0.10 


0.34 


0.19 


0.29 


0.08 


0.31 


0.08 


0.41 


0.11 


0.36 


0.08 


0.39 


0.07 


0.54 


0.22 






0.54 


0.22 






0.21 


0.07 


0.34 


0.07 


0.23 


0.10 


0.31 


0.08 


0.20 


0.08 


0.15 


0.09 


0.35 


0.13 


0.13 


0.07 


0.07 


0.07 


0.16 


0.08 


0.16 


0.09 


0.12 


0.08 


0.35 


0.06 


0.41 


0.06 


0.30 


0.06 


0.33 


0.06 


0.24 


0.22 






0.24 


0.22 






0.25 


0.22 






0.25 


0.22 






0.31 


0.18 






0.31 


0.18 


0.32 


0.19 


0.22 


0.07 


0.21 


0.06 


0.29 


0.11 


0.33 


0.06 


0.25 


0.08 


0.29 


0.06 


0.15 


0.07 


0.41 


0.07 


0.31 


0.06 


0.30 


0.09 


0.40 


0.06 


0.25 


0.06 


0.39 


0.07 


0.29 


0.06 


0.42 


0.10 


0.34 


0.07 


0.14 


0.18 






0.14 


0.18 






0.20 


0.07 


0.27 


0.07 


0.15 


0.08 


0.26 


0.07 


0.33 


0.07 


0.28 


0.08 


0.44 


0.09 


0.31 


0.07 


0.13 


0.11 


0.21 


0.06 


0.15 


0.12 


0.21 


0.06 


0.33 


0.07 


0.19 


0.06 


0.44 


0.12 


0.28 


0.07 


0.35 


0.05 


0.41 


0.07 


0.19 


0.13 


0.42 


0.07 


0.26 


0.09 


0.22 


0.07 


0.31 


0.11 


0.22 


0.07 


0.23 


0.07 


0.16 


0.06 


0.19 


0.08 


0.21 


0.06 


0.20 


0.10 


0.31 


0.10 


0.20 


0.15 


0.29 


0.14 


0.21 


0.15 


0.39 


0.07 


0.31 


0.21 


0.29 


0.08 


0.34 


0.08 


0.30 


0.07 


0.30 


0.11 


0.28 


0.07 


0.26 


0.20 


0.42 


0.06 


0.25 


0.17 


0.41 


0.06 


0.23 


0.15 


0.28 


0.07 


0.11 


0.21 


0.35 


0.07 


0.27 


0.11 


0.35 


0.07 


0.36 


0.21 


0.37 


0.07 


0.41 


0.19 


0.25 


0.08 


0.28 


0.22 


0.27 


0.08 


0.25 


0.18 


0.27 


0.08 


0.21 


0.22 


0.19 


0.08 


0.38 


0.21 


0.27 


0.12 


0.28 


0.22 


0.36 


0.11 


0.23 


0.08 


0.27 


0.06 


0.28 


0.21 


0.29 


0.06 


0.38 


0.12 


0.36 


0.10 


0.44 


0.15 


0.35 


0.09 


0.07 


0.22 


0.30 


0.13 


0.27 


0.14 


0.30 


0.12 



a In/out split at d ~ 12.1 kpc 
b In/out split at d ~ 14.1 kpc 
c In/out split at d ~ 16.8 kpc 
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Table 3 

Mean [Fc/H] in the Smooth Component of the Halo 



Volume 


Mean [Fc/H] 


err 


Mean [Fe/H] Above |6| = 30° 


err 


10.0 kpc < d < 12.1 kpc 


-1.77 


0.12 


-1.78 


0.12 


10.0 kpc < d < 14.1 kpc 


-1.78 


0.12 


-1.79 


0.11 


10.0 kpc < d < 16.8 kpc 


-1.72 


0.12 


-1.77 


0.12 


10.0 kpc < d < 17.5 kpc 


-1.65 


0.13 


-1.75 


0.12 


12.1 kpc < d < 17.5 kpc 


-1.58 


0.13 


-1.69 


0.12 


14.1 kpc < d < 17.5 kpc 


-1.55 


0.13 


-1.63 


0.12 


16.8 kpc < d < 17.5 kpc 


-1.45 


0.13 


-1.53 


0.13 
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Table 4 

Spatial Autocorrelation in [Fe/H] in the Smooth Component of the Halo 



Volume 


I 


E(J) 


Var(/) 


p- value 


10.0 kpc < d < 12.1 kpc 


-0.012 


-0.012 


0.019 


l.Oc+00 


10.0 kpc < d< 14.1 kpc 


0.018 


-0.010 


0.016 


8.2e-02 


10.0 kpc < d < 16.8 kpc 


0.049 


-0.009 


0.015 


1.2e-04 


12.1 kpc < d < 17.5 kpc 


0.110 


-0.009 


0.015 


4.4e-16 


14.1 kpc < d < 17.5 kpc 


0.113 


-0.009 


0.016 


3.1c-15 


16.8 kpc < d < 17.5 kpc 


0.077 


-0.010 


0.018 


1.3e-06 



Note. — The five columns are: range in distance from the Galactic center, Moran's I statistic indicating the degree of spatial 
autocorrelation in the data in Tabic [l] the expected value of Moran's / statistic under the null hypothesis that there is no spatial 
autocorrelation, the variance of the expected value under the null hypothesis, and the probability that the null hypothesis applies to the 
data in Table [T] 
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Table 5 

Spatial Autocorrelation in [Fe/H] Above \b\ = 30° in the Smooth Component 

of the Halo 



Volume 




E(J) 


Var(7) 


p- value 


10.0 kpc < d < 12.1 kpc 


-0.033 


-0.013 


0.021 


3.5e-01 


10.0 kpc < d < 14.1 kpc 


-0.018 


-0.011 


0.018 


6.7e-01 


10.0 kpc < d < 16.8 kpc 


0.026 


-0.010 


0.017 


3.5e-02 


12.1 kpc < d < 17.5 kpc 


0.070 


-0.010 


0.018 


4.7e-06 


14.1 kpc < d < 17.5 kpc 


0.088 


-0.011 


0.019 


l.lc-07 


16.8 kpc < d < 17.5 kpc 


0.067 


-0.012 


0.022 


2.1c-04 
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Tabl e 6 

Spatial Autocorrelation in [Fc/H] in Bullock & Johnston (2005) Halos 

Halo Volume / Var(J) p-value 



1 


10.0 kpc < d < 12.1 kpc 


0.006 


-0.009 


0.015 


3.3c-01 




10.0 kpc < d < 14.1 kpc 


0.001 


-0.008 


0.014 


5.1c-01 




10.0 kpc < d < 16.8 kpc 


-0.016 


-0.008 


0.015 


5.9e-01 




12.1 kpc < d < 20.0 kpc 


-0.009 


-0.008 


0.015 


9.6e-01 




14.1 kpc < d < 20.0 kpc 


0.044 


-0.008 


0.015 


4.3e-04 




16.8 kpc < d < 20.0 kpc 


-0.023 


-0.038 


0.048 


7.4e-01 



10.0 kpc < d < 12.1 kpc 0.011 -0.009 0.016 2.1e-01 

10.0 kpc < d < 14.1 kpc -0.005 -0.009 0.015 8.3e-01 

10.0 kpc < d < 16.8 kpc 0.008 -0.008 0.015 2.8e-01 

12.1 kpc < d < 20.0 kpc 0.018 -0.008 0.015 7.8e-02 
14.1 kpc < d < 20.0 kpc 0.015 -0.008 0.015 l.lc-01 
16.8 kpc < d < 20.0 kpc 0.009 -0.030 0.059 5.0c-01 
10.0 kpc < d < 12.1 kpc -0.016 -0.009 0.014 6.1c-01 
10.0 kpc < d < 14.1 kpc -0.005 -0.008 0.013 7.8e-01 

10.0 kpc < d < 16.8 kpc 0.002 -0.008 0.014 4.7e-01 

12.1 kpc < d < 20.0 kpc 0.018 -0.008 0.015 7.6e-02 
14.1 kpc < d < 20.0 kpc 0.024 -0.009 0.015 3.3e-02 
16.8 kpc < d < 20.0 kpc -0.047 -0.037 0.066 8.8e-01 
10.0 kpc < d < 12.1 kpc 0.018 -0.009 0.016 8.6e-02 
10.0 kpc < d < 14.1 kpc 0.002 -0.008 0.015 4.7e-01 

10.0 kpc < d < 16.8 kpc 0.012 -0.008 0.015 1.7e-01 

12.1 kpc < d < 20.0 kpc -0.001 -0.008 0.015 6.0e-01 
14.1 kpc < d < 20.0 kpc 0.004 -0.008 0.015 4.0e-01 
16.8 kpc < d < 20.0 kpc 0.264 -0.023 0.037 1.2e-14 
10.0 kpc < d < 12.1 kpc 0.005 -0.009 0.016 3.9e-01 
10.0 kpc < d < 14.1 kpc -0.002 -0.009 0.015 6.9e-01 

10.0 kpc < d < 16.8 kpc -0.005 -0.009 0.015 8.0e-01 

12.1 kpc < d < 20.0 kpc -0.003 -0.009 0.015 7.0e-01 
14.1 kpc < d < 20.0 kpc -0.009 -0.009 0.015 9.9e-01 
16.8 kpc < d < 20.0 kpc -0.024 -0.028 0.050 9.3e-01 
10.0 kpc < d < 12.1 kpc 0.037 -0.009 0.016 3.9e-03 
10.0 kpc < d < 14.1 kpc 0.025 -0.008 0.015 2.6e-02 

10.0 kpc < d < 16.8 kpc 0.037 -0.008 0.015 2.5e-03 

12.1 kpc < d < 20.0 kpc 0.006 -0.008 0.015 3.3e-01 
14.1 kpc < d < 20.0 kpc 0.000 -0.008 0.014 5.4e-01 
16.8 kpc < d < 20.0 kpc -0.052 -0.023 0.038 4.5e-01 
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10.0 kpc < d < 12.1 kpc 


-0.012 


-0.009 


0.016 


8.8c-01 




10.0 kpc < d < 14.1 kpc 


0.021 


-0.009 


0.015 


5.1c-02 




10.0 kpc < d < 16.8 kpc 


0.037 


-0.008 


0.015 


2.7e-03 




12.1 kpc < d < 20.0 kpc 


0.016 


-0.008 


0.015 


1.0e-01 




14.1 kpc < d < 20.0 kpc 


0.000 


-0.009 


0.015 


5.6e-01 




16.8 kpc < d < 20.0 kpc 


0.067 


-0.036 


0.064 


l.lc-01 



10.0 kpc < d < 12.1 kpc -0.007 -0.009 0.016 9.1e-01 

10.0 kpc < d < 14.1 kpc -0.015 -0.008 0.015 6.7e-01 

10.0 kpc < d < 16.8 kpc -0.004 -0.008 0.015 7.5e-01 

12.1 kpc < d < 20.0 kpc 0.022 -0.008 0.015 3.8e-02 
14.1 kpc < d < 20.0 kpc 0.026 -0.008 0.015 1.9e-02 
16.8 kpc < d < 20.0 kpc 0.042 -0.018 0.032 6.5e-02 
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10.0 kpc < d < 12.1 kpc 


-0.022 


-0.009 


0.016 


4.2e-01 




10.0 kpc < d < 14.1 kpc 


-0.035 


-0.008 


0.015 


7.6e-02 




10.0 kpc < d < 16.8 kpc 


-0.006 


-0.008 


0.015 


8.6e-01 




12.1 kpc < d < 20.0 kpc 


0.019 


-0.008 


0.014 


5.6e-02 




14.1 kpc < d < 20.0 kpc 


0.021 


-0.008 


0.015 


4.5e-02 




16.8 kpc < d < 20.0 kpc 


-0.034 


-0.029 


0.034 


8.6e-01 
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10.0 kpc < d < 12.1 kpc 


0.002 


-0.009 


0.016 


5.0c-01 




10.0 kpc < d < 14.1 kpc 


0.021 


-0.008 


0.014 


4.1c-02 




10.0 kpc < d < 16.8 kpc 


-0.002 


-0.008 


0.014 


6.6e-01 




12.1 kpc < d < 20.0 kpc 


-0.007 


-0.008 


0.015 


9.5e-01 




14.1 kpc < d < 20.0 kpc 


-0.023 


-0.008 


0.015 


3.3e-01 




16.8 kpc < d < 20.0 kpc 


-0.055 


-0.028 


0.037 


4.7c-01 
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10.0 kpc < d < 12.1 kpc 


0.041 


-0.014 


0.024 


2.2e-02 




10.0 kpc < d < 14.1 kpc 


0.036 


-0.011 


0.019 


1.2e-02 




10.0 kpc < d < 16.8 kpc 


0.027 


-0.009 


0.016 


2.5e-02 




12.1 kpc < d < 20.0 kpc 


0.036 


-0.008 


0.015 


2.6e-03 




14.1 kpc < d < 20.0 kpc 


0.057 


-0.009 


0.016 


2.4e-05 




16.8 kpc < d < 20.0 kpc 


0.049 


-0.010 


0.018 


9.7e-04 



THE FORMATION OF THE MILKY WAY THROUGH COLD HALO SUBSTRUCTURE. 

Table 7 

Spatial Autocorrelation in [Fe/H] in Via Lactea 2 



Volume / Var(7) p-value 

10.0 kpc < d < 12.1 kpc 0.025 -0.019 0.033 1.8c-01 

10.0 kpc < d < 14.1 kpc 0.067 -0.015 0.024 7.5e-04 

10.0 kpc < d < 16.8 kpc 0.003 -0.012 0.020 4.4e-01 

12.1 kpc < d < 20.0 kpc 0.079 -0.011 0.019 1.3e-06 
14.1 kpc < d < 20.0 kpc 0.078 -0.012 0.020 8.9e-06 
16.8 kpc < d < 20.0 kpc 0.043 -0.015 0.025 1.9e-02 



